Submitted to Ap J L 



Gas Density, Stability and Starbursts Near the 
Inner Lindblad Resonance of the Milky Way 

Antony A. Stark, Christopher L. Martin, Wilfred M. Walsh, Kecheng Xiao, Adair P. Lane 

Harvard- Smithsonian Center for Astrophysics, 60 Garden St., Cambridge MA 02138 

aasOcf a. harvard. edu 

and 

Christopher K. Walker 

Steward Observatory, University of Arizona, Tucson, AZ 85721 
cwalkerOas . arizona . edu 

ABSTRACT 

A key project of the Antarctic Submillimeter Telescope and Remote Obser- 
vatory reported by Martin et al. (2004) is the mapping of CO J = 4 — >^ 3 and 
J = 7 — >■ 6 emission from the inner Milky Way, allowing determination of gas 
density and temperature. Galactic center gas that Binney et al. (1991) identify as 
being on X2 orbits has a density near 10^'^ cm~^, which renders it only marginally 
stable against gravitational coagulation into a few Giant Molecular Clouds, as 
discussed by Elmegreen (1994). This suggests a relaxation oscillator mechanism 
for starbursts in the Milky Way, where inflowing gas accumulates in a ring at 150 
pc radius for approximately 20 million years, until the critical density is reached, 
and the resulting instability leads to the sudden formation of giant clouds and 
the deposition of 4 x lO^M© of gas onto the Galactic center. 



Subject headings: Galaxy: structure — ISM: clouds — ISM: molecules — galaxies: 
st arburst — st ars : formation 
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1. Introduction 

Dynamics of gas in the inner few kiloparsecs of the Milky Way are dominated by the 
non-axisymmetric gravitational potential of the central bar. The properties of this bar are 
now well known, and there is good agreement between observations of gas motion and model 
fits to the potential (Jenkins & Binney 1994; Gerhard 1999; Hafner et al. 2000; Bissantz 
et al. 2003). As suggested by Binney et al. (1991), the gas tends to be found on famihes of 
closed orbits which are not self-intersecting. All non-closed orbits and some closed orbits are 
self-intersecting. Gas on such orbits will shock and lose energy where the gas streamlines 
intersect, and the gas will then move inwards to a lower energy orbit. If the gas can find its 
way onto a family of non-self-intersecting closed orbits, the energy dissipation slows and the 
timescale for orbital changes lengthens out. Contopoulos & Mertzanides (1977) described 
two families of closed orbits in barred galaxies: the "xi", which are elongated along the 
bar and found outside the inner Lindblad resonance (ILR); and the "X2", which are more 
round and can be found near the ILR and inside it. The ILR is located where the epicyclic 
frequency of a particle orbiting in the Galactic potential resonates with the pattern speed of 
the bar. This occurs at a radius of approximately 450 pc from the center of the Milky Way 
(Bissantz et al. 2003). 

Gas which is several kiloparsecs away from the Galactic center tends to be driven inwards 
until it reaches a region within a few hundred parsecs of the ILR, because the interaction of 
the bar potential with the gas exerts a negative torque, resulting in loss of angular momentum 
by the gas (Lynden-Bell & Kalnajs 1972; Athanassoula 1988; Jenkins & Binney 1994). Near 
the ILR this effect disappears, because the net torque there is small or zero, and inwards 
of the ILR it may even reverse and become positive, so that gas inside the ILR could be 
driven outwards (Combes 1988). Gas therefore accumulates in a ring near the ILR. Unlike 
the gas further out, the dynamics of this gas depends critically on its self gravity (Elmegreen 
1994; Jenkins & Binney 1994), and therefore on its thermodynamic properties, density in 
particular. 

The thermodynamic properties of the gas can be determined by millimeter- and sub- 
millimeter-wave spectral line observations. The distribution of molecular gas near the ILR 
is known from extensive surveys in CO and ^^CO J = 1 — > and J = 2 — >■ 1 (Bally 
et al. 1988; Bitran et al. 1997; Oka et al. 1998); these spectral lines show the presence 
of molecular gas. These hnes alone do not, however, determine its density or excitation 
temperature. Observations of the mid- J lines of CO provide the missing information. Since 
the low- J states of CO are in local thermodynamic equilibrium (LTE) in almost all molecular 
gas (Goldreich & Kwan 1974), measurements of mid- J states are critical to achieving a 
solution of the radiative transfer by breaking the degeneracy between beam filling factor and 
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excitation temperature. A new survey (Martin et al. 2004) by the Antarctic Submillimeter 

Telescope and Remote Observatory (AST/RO, Stark et al. 2001) adds the J = 4 — >■ 3 
(461 GHz) and J = 7 ^ 6 (807 GHz) rotational lines of CO to existing lower- frequency 
data (Bally et al. 1988). These data are available on the AST/RO website^ for general use. 
These measurements have recently been modeled using the large velocity gradient (LVG) 
approximation to determine the gas density and temperature. 

In this Letter, we discuss the implications of the Martin et al. (2004) density estimates 
for Galactic center gas. We apply our new data, specific to the Milky Way, to the general 
analysis of stability of dense gas near ILR regions in galaxies by Elmegreen (1994). We find 
that the gas near 150 pc radius is marginally unstable. This suggests that in the past there 
has been a period of stability and gas build-up. In the future, the instability will create a few 
giant clouds, resulting in a starburst and the deposition of tens of millions of solar masses 
of material on the Galactic center. This process repeats with a cycle time determined by 
the rate at which gas precipitates on the Galactic center region from outside, resulting in 
starbursts at intervals of approximately 20 million years. 



2. Gas Density 

Figure 1 is a representation of the average density of the molecular gas layer in the 
vicinity of the ILR. Martin et al. (2004) used an LVG model on their survey data to estimate 
the density and temperature at each point in (£, b, v). Density values in Figure 1 are 
calculated by averaging the logarithm of the density values from Martin et al. (2004) at 
each value of {i.v), as b varies from —0.3 to 0.2, excluding points where the LVG fit is 
uncertain. The excluded points are those where the average ^^CO J = 1 ^ emission in 
Bally et al. (1988) is less than 0.2K, or where the average density is below 10^'^ cm~^, which 
is an approximate lower limit to the validity of the Martin et al. (2004) LVG model. The 
pixels at some values of {£, v) in Figure 1 appear white, and this can be for several reasons: 
(1) there are no data at that position since the survey is limited in v by the bandwidth 
of the spectrometer; (2) there is no significant CO emission; (3) the density is below the 
threshold of validity of the model; (4) the variance in the average over b exceeds (10^ cm~^)^. 
The model has converged to a consistent value at the points where a density is displayed. 
Even those values could be spurious if the assumptions used in the LVG approximation do 
not apply. This certainly occurs where there is foreground absorption far from the emitting 
region, for example between —60 and +25kms^^, where spiral arms in the Galactic disk are 
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projected onto the Galactic center emission. This region is hatched in Figure 1. 

Superposed on the data are some closed orbits from the model of Bissantz et al. (2003). 
The cusped xi orbit and two self-intersecting xi orbits are shown in magenta. The cusped 
xi orbit is the innermost xi orbit which does not self- intersect, with an apogee 1300 pc from 
the Galactic center. The two self-intersecting xi orbits are interior to this, with apogees 
near 1000 pc. The X2 orbits begin inside the ILR, and are drawn in brown. The outermost 
X2 has an apogee of 180 pc. The X2 orbits extend all the way into the Galactic center; the 
innermost shown here is almost circular and has an apogee of 23 pc. These same orbits are 
plotted differently in Figures 10 and 11 of Bissantz et al. (2003). 

We see that in some places along the cusped xi orbit, n{B.2) ~ 10^ to 10"^ cm~^, although 
the gas is clumpy and the average density on this orbit is less. Most places along the outer X2 
orbit have n(H2) ~ 2 x 10^ to 6 x 10'^ cm~^, with few values outside this range. The densities 
on the X2 orbits which are at negative £ and negative v are unreliable, because these parts of 
the orbit lie about 1 kpc behind some of the outer xi orbits (not shown in Figure 1) which 
have the same velocity but lower excitation. Most of the region between the cusped Xi orbit 
and the outermost X2 orbit has lower density and little molecular emission. 



3. Gas Stability 

Given a measure of the gas density on orbits near the ILR, it is possible to determine 
whether or not that gas is stable — if the gas is sufficiently dense, its self-gravity will overcome 
the tidal shear of the Galaxy's gravitational potential and it will agglomerate into clouds. 
The dynamics of this situation has been analyzed by Elmegreen (1994), who hnearized the 
hydromagnetic force and continuity equations in a rotating frame and thereby derived a dis- 
persion relation for the growth of instabilities in a gas ring near the ILR. He determined that 
the growth rate is a function of gas density and pressure, the Galactic potential, magnetic 
field strength, and the rate at which material accretes from larger radii. This relation is 
expressed in equation 11 of Elmegreen (1994): 



where ujr is the growth rate of instability, k is its wavenumber, k is the epicyclic frequency 
in the Galactic potential, va is the Alfven velocity corresponding to the azimuthal magnetic 
field, Qa is the relative accretion rate, and ujg is a frequency related to the acceleration of 
self-gravity. Values of the quantities in Equation 1 are estimated in Table 1. 



The relative accretion rate, Qa, is a measure of gas inflow from larger radii. It is 
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determined by processes outside the ILR region: the amount of gas in the outer regions 
of the bar, and the torques exerted on that gas by the bar. At minimum, the evolved 
stars in the bulge will eject matter into the interstellar medium of the outer bar at a rate 
~ 0.2MQyr~^ (Jungwiert et al. 2001). At maximum, the Galactic center region could ingest 
an entire gas-rich dwarf companion at a rate of 100 to 1000 MQyr^^. Accretion rates higher 
than this would disrupt the inner Galaxy, so we can conclude they have not occurred in the 
Milky Way (Heller & Shlosman 1994; Bournaud & Combes 2002). Combes (2004) suggests 
an accretion rate of lOMoyr"^, which is approximately the rate of accretion of intergalactic 
material onto the Galaxy. This rate, accreting onto inner gas disk of ~ 2 x 10^ M©, yields 
Qa ~ SOGyr-^ 

There is as yet no measurement of the strength of the magnetic field in the dense Galac- 
tic center gas, although the existence of a magnetic field perpendicular to the disk near 
the Galactic center is demonstrated by the non-thermal radio filaments (Yusef-Zadeh et al. 
1984), and the existence of a magnetic field in the plane of the disk is demonstrated by 
submillimctcr-wave polarimetry (Novak et al. 2003). Chuss et al. (2003) have estimated the 
field strength to be about a milliGauss, based on the field morphology. The perpendicular 
field near the center, they argue, is the result of a process where the field in the disk is 
amplified until the magnetic field pressure begins to dominate the dynamics and the field 
decouples from the gas in the regions of relative lower density. This implies that the magnetic 
field energy in the regions of relative higher density self-regulates to approximate equiparti- 
tion with the internal kinetic energy of the gas, making the Alfven velocity, va, comparable 
to the internal turbulent velocity of the molecular gas. Accepting this argument leads to a 
value of Va in Equation 1 which is significant but not dominant. In other barred galaxies, 
the observed field strength averaged over the central few kiloparsecs is typically 10 fiG (Beck 
et al. 2002). Since the dense gas has a filling factor of ~ 10~^, this is consistent with a 
milliGauss field in the dense gas. 

The epicyclic frequency, k, is more certain. It depends only on the rotation curve and 
its derivative, which are known within ±20%. In Table 1 we adopt values taken from Figure 
1 of Bissantz et al. (2003). 

Elmegreen (1991) approximates the equation of state for the molecular gas as AP f« 
7efrc^Ap, where the complexities of heating and cooling of the molecular gas are subsumed 
in an effective ratio of specific heats 'jcS ~ 0.3 to 2, and an effective sound speed c ~ 10 to 
40kms~^. As pointed out in Elmegreen (1994), the quantity 7cfrc^ for the actual gas under 
analysis can be estimated by considering the equilibrium between pressure and self-gravity 
perpendicular to the plane: h ~ {2 / e)'^^^ ck~^ , where the scaleheight, h, is determined 
from observations. Applied to the values in Table 1, this procedure gives 7eff ^ 1 if c ~ ~ 
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25kms-^ 

The gravitational frequency, lug: as defined by Elmegreen (1994) is a measure of the 
relative importance of self-gravity compared to pressure in the gas over a scalelength given 
by wavenumber k: 

a;© ~ 2 Gi^e In (^^^ - P -f.sc' (2) 

where fi = pnh^ is the mass per unit length in the ring. The uoq has a maximum when 
^"max ~ 2/i^^exp[— 0.5(1 + 7cfrC^/G/i)], and the corresponding gravitational frequency, ujq ~ 
G/^^max corresponds to the fastest growing unstable mode. 

We are now in a position to evaluate u^. the growth rate of the fastest unstable mode, 
by substituting /Cmax for k. Equation 1 is quartic in uj^ with one positive real solution. If 
ila <C and k^^VA ^ as we see is the case from Table 1, then 

cUr <^ ojg if ojg ^ and 

UJr ~ OJg if OJg ^ 1^ ■ 

In the Milky Way at the present time, the accretion term Via is not dynamically significant to 
the instability (although it does determine the time between instabilities), and the magnetic 
field term va has only marginal significance, since its contribution to the dynamics is in 
approximate equipartition with the other terms. The growth rate of the instability is small 
as long &s ujg ^ but becomes significant when ujg ~ The dominant dynamical effect 
is the competition between the tidal shear of the Galaxy and the tendency for the gas ring 
to clump with wavelength Amax = 27r/A;max- The criterion for significant instability, ujg > k, 
can be expressed as a threshold in gas density: 

(3) 

which is essentially a density criterion for the formation of molecular clouds in the presence 
of Galactic tidal forces (Stark et al. 1989). The region that begins to contract will initially 
have a mass approximately equal to Umax^max = 4 x 10^ M0. As the instabihty proceeds, the 
dynamics will become non- linear, and move beyond the valid regime of Equation 1. This 
merits further study. 

4. Starbursts in the Milky Way 

We see from Figure 1 that gas on the inner xi orbits exceeds the critical density threshold 
in places, and that the gas is clumped. The gas on the outer X2 orbits is more smooth, but 
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is at the threshold of significant instabihty. The timescale for clumping to begin is short: 
1/ojr ~ 1 Myr . Since Amax ~ r, we expect the instability to result in only a few large clouds, 
containing many millions of solar masses. The Giant Molecular Cloud surrounding Sgr B2, 
which can be seen as a relatively minor density enhancement in Figure 1 at £ = 0.65°, 
V = 60kms~^, is such a cloud. Because of their large mass, these clouds arc subject to 
dynamical friction with the background of stars in the Galactic bulge, and will spiral into 
the center within 1 Gyr (Stark et al. 1991). Dynamical friction can overcome the forces 
tending to maintain the gas in the ILR region, and allow the gas to continue inwards toward 
the galactic center, but only if the gas is organized into sufficiently large self-gravitating 
clouds. This suggests a relaxation oscillator mechanism for quasi-periodic starbursts in the 
center of the Milky Way. At first, the Galactic center region is relatively clear of gas, like 
the central region of M31 is now. Gas precipitates into the region of the bar, either as 
mass loss from evolved stars or as the result of cannibalism of smaller galaxies. The bar 
dynamics drive this gas toward the ILR, where it will tend to accumulate in a ring as long 
as n(H2) ^ ncrit- The ring becomes more and more dense as gas continues to precipitate 
from larger radii, and eventually the threshold is reached. A few giant clouds will form on 
a relatively short timescale, creating a star burst, and the giant clouds will be swept toward 
the center by dynamical friction, restoring a condition of relatively low density. The cycle 
will repeat on timescales of ~ 20 Myr, but of coTirsc this timescale is highly variable 
and can be dramatically shortened by events which precipitate a large amount of gas. 

Support was provided by NSF grant OPP-0126090. We thank our AST/RO colleagues 
Richard Chamberlin, Jacob Kooi, and Gregory Wright for extensive help with the instru- 
mentation. 
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Table 1. Quantities relating to the stability of gas near the ILR. 



Quantity 




Value at r = 150 pc 


Value at r = 450 pc 


Accuracy 


Reference 


iVitix. inoicculcir miiriDGr clGiisity 


'tmax {1^2 } 




10^ cm~^ 




ividJ Lin et ai. ^zuu^^ 


Al(ix. ^tis density 


Pmax — mji 'lmaxl,ri2^ 










Rsltitiv© cLCcrstioii rcitc 


o 




ou i-jyi 




ooniucs izuu-i) 


Azimuthal magnetic field 


Be 


1 mO 


1 mG 


C 


Chuss et al. (2003) 


Az. Alfven velocity 


VA = Be{47rpmax)"^''^ 


26kms~l 


20kms~l 


c 




Galactic angular velocity 


n 


620 Gyr-1 


314Gyr-l 


A 


Bissantz et al. (2003) 


Derivative of Q 


dn/dr- 


-2500Gyr-ikpc-i 


-490 Gyr-^kpc-l 


A 


Bissantz et al. (2003) 


Epicyclic frequency 


K = (4^2 + 2rQdn/dr-)V2 


1036 Gyr-i 


506 Gyr-i 


A 




Molecular half-scaleheight 


h 


21 pc 


56 pc 


A 


Bally et al. (1988) 


Ratio of specific heats 


7effc2 = AP/Ap « (/i«e)V4 


874 km^ s"2 


1483 km^s-^ 


B 


Elmegreen (1991) 


Max. mass per length of ring 


/imax — T^Pmaxh^ 


2.5 X lO^MQkpc-i 


2.8 X lO^Mokpc-i 


B 


Elmegreen (1994) 


Wavenumber of max. growth 


/?max « 2/l-lexp[-0.5(l + 7effC^/<5/imax)] 


38kpc-i 


20kpc-l 


B 


Elmegreen (1994) 


Wavelength of max. growth 


-^max — 27r//Cniax 


165 pc 


314 pc 


B 


Elmegreen (1994) 


Gravitational frequency at fcmax 


= fcmax(G/imax)"'"''^ 


1270 Gyr-i 


2286 Gyr^l 


B 


Elmegreen (1994) 


Instability growth rate 


OJr 


1013 Gyr-i 


2206Gyr-i 


B 





Note. — The Sun's distance to the Galactic center is taken to be Rq = 8 kpc. The effective sound speed in the interstellar gas is c. Galactocentric radius is r. The 
column labeled "Accuracy" indicates the approximate errors — A: ±20%; B; ±50%; C: order-of-magnitude. 
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Fig. 1. — Density on xi and X2 orbits in the inner Milky Way. The rainbow scale at right 
shows the average density of molecular gas from an LVG model using AST/RO survey data 
(Martin et al. 2004). White pixels are points where there are no data, or where the LVG 
model does not converge to consistent values. The hatched area at low velocity shows definite 
foreground absorption which invalidates the assumptions of the LVG model. Superposed on 
the density data are some xi (magenta) and X2 (brown) orbits from Bissantz et al. (2003). 



